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Abstract 



Assuming the existence of stationary metastable states (SMS's), and 
using general thermodynamic arguments, we prove that a positive Kauz- 
mann temperature exists for SMS's provided the ideal glass energy is 
higher than the crystal energy at absolute zero. We confirm the general 
predictions by two exact calculations, one of which is not mean-field. 

Glass transition is a ubiquitous phcnomcnon^'^ which is believed to occur 
in a supercooled liquid (SCL), and has been investigated for at least over eight 
decades since the earliest works of Nernst and Simon'^. Despite this, a complete 
understanding of the transition and related issues is still far from complete, al- 
though major progress has been made recently.^"^ Theoretical and experimental 
investigations invariably require applying (time-independent) thermodynamics 
to SCL to extract quantities like the entropy. This presumes, as we do here, 
that there exists a stationary limit (SMS) of the metastable state (SCL) under 
infinitely slow cooling in which the crystal is forbidden from nucleating.^ We 
further assume that SCL can be defined all the way down to absolute zero, 
which may not be the case always.^ 

We focus on the partition function (PF) determined by the configurational, 
i.e., the potential energy and the related entropy in this work due the cen- 
tral importance of the configurational entropy, a classical statistical mechanics 
(CSM) concept,^ for the ideal glass transition. Experimentally, the SCL con- 
figurational entropy exhibits a rapid drop near a temperature ~ two-thirds of 
the melting temperature Tm,^'^ and its smooth low-temperature extrapolation 
gives rise to a negative entropy (the entropy crisis commonly known as the 
Kauzmann paradox)^ below a non-zero Kauzmann temperature Tk. There are 
several computational^ and theoretical^'^^^*^ results clearly demonstrating the ex- 
istence of the entropy crisis. A positive Tk is also conceptually necessary for the 
observed Vogel-Fulcher behavior in fragile systems. However, the situation has 
become very confusing, as there have appeared several recent arguments"'^^"^^ 



against its existence. It is, therefore, extremely important to clarify the is- 
sue, which we do by proving in two different ways the existence of a non-zero 
Tk using general but rigorous thermodynamic arguments, valid for classical or 
quantum systems. The proof assumes SMS- existence, so that thermodynamics 
can be applied. Its stationary nature allows us to investigate SMS using the PF 
formalism in (1) below. We verify our conclusions by two exact calculations, 
one of which is not a mean-field calculation. 

Consider a system composed of N particles for which E, V and S, which are 
homogeneous functions of order 1 in N, denote the energy, the volume, and the 
entropy. The canonical PF for fixed V and N is given by 

Z = TrW{E) exp{-l3E), (1) 

where Tr stands for summation over the energy (eigenvalues) E > Eq up to 
its maximum, /? = 1/T, T being the system temperature in the units of the 
Boltzmann constant fcs, and W{E) = exp(S') > 1 is defined suiatbly to rep- 
resent the number of distinct microstates for a given E. The lowest allowed 
energy Eq [W{Eo) ^ 0] certainly represents the Hclmholtz free energy and the 
energy of the perfect crystal (CR) at T = 0. In equilibrium, Z must be max- 
imized in the thermodynamic limit N — > oo, keeping E/N,V/N fixed. The 
entropy^ S{E) = In W{E) must be non-negative so that the state is observable 
in Nature.^"' This observability condition on S cannot be violated, even for a 
SMS if it is observable. Therefore, its violation (SsMsiT) < for T < Tk > 0)" 
implies that SMS cannot exist in Nature when the violation occurs, in which 
case the SMS associated with the SCL must be replaced by a new state, com- 
monly known as the ideal glass state below Tk, whose energy at T < Tk > 
is Ek > Eq}'"^ The transition between the two states is called the ideal glass 
transition. 

The order parameter p (see also below) is traditionally used to distinguish 

various phases of the system, with p = representing the disordered phase, 
the equilibrium liquid (EL), containing the number of microstates W^is consis- 
tent with p = 0. The ordered phase CR has p ^ 0, and contains the num- 
ber of microstates Word consistent with p 7^ 0. We partition W(B) using p : 
W{B) = VFord(-B) -|- Wdis(-B). (Presence of other configurations will not affect 
our argument.) We introduce corresponding PF Z„ (a = dis or ord) following 
(1) with W replaced by Wa- Above Tm, Zdis >^ord, and below Tm, Zdis<^ord 
. The PF along with the its maximization requirement allows us to describe 
CR and EL collectively. We describe SMS as the continuation of Z^ig below 
Tm by abandoning this requirement. This necessarily yields p = also for 
SCL. While Wa{E) > 1 for physically observed states, there is no reason for it 
to hold everywhere for SMS states obtained by the above continuation; states 
with Wa{E) < 1 are not realizable. -'^^ Nevertheless, as long as Wa{E) > 0, Z^ 
is a sum of positive terms, so that fla = InZ^, (i) satisfies proper convexity 
properties'^, and (ii) is dominated by the maximum value of Sa — PE in Z^. 
Expressing 

Z« = W„(£;o)e-^^°(l + Z„), (2) 
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where = Tr \y/a{E)IWa{E^^)Y'^^^^'^"^ ■, with Tr restricted to all E > Eq, 
we find that Tn^ = TSa{Eo) - Eq + rin(l + Z^). For all finite N,Za-^0 
as T — > 0, provided Wa{Eo) > 0. The last condition means that TSa — * 
as T ^ 0: hence. F,, = -TO,, Eq for CR and SMS at T = : Both have 
identical free energies at absolute zero, so that F^/ Eq = 1. Now, we can take the 
limit N ^ oo without altering the conclusion Fa/Eo = 1. However, it is known 
experimentally that (1) the heat capacity of the glassy state can be substantially 
different from that of the corresponding crystal at the same temperature,^^ and 
(2) the energy of this state Ek extrapolated to absolute zero is higher than that 
of the CR energy Eq}''^ Consequently, the stationary limit of the glassy state is 
different not only from the CR, but also from the SMS near absolute zero. Since 
SMS has a positve heat capacity, it achieves its energy Ek > Eq at a positive 
temperature Tk, where it merges with the ideal glass. This is our first rigorous 
proof that an ideal glass transition must occur at a positive temperature if we 
identify the stationary limit of the glassy state with the ideal glass noted above. 

We now proceed to give an alternative rigorous proof independent of Fa/ Eq = 
1 at T = derived in the canonical ensemble. We consider the microcanonical 
ensemble in which the entropy functions Sa for fixed V,N, are shown schemati- 
cally as functions of E by OAB for the CR (OA) and its SMS extension (AB) for 
the superheated CR, and DACK for the EL (DA), and its SCL extension (ACK) 
in Fig. 1. The corresponding Helmholtz free energies F^ vs T in the canonical 
ensemble are shown schematically by BO and DACK in the inset. Since a SMS 
is not an equilibrium state, its entropy at some E cannot exceed the entropy 
of the corresponding equilibrium state at the same E. This explains the form 
of the entropy in Fig. 1, which is also supported by all known observations, 
exact®'^'^'^'^ and numerical calculations,^ and from the arguments and the calcu- 
lations presented below. We also show the entropy function for a non-stationary 
metastablc state by OF in Fig. 1, assuming that crystallization is forbidden. 
The slope of HH gives I/Tm and of 00' gives the inverse temperature (see O 
in the inset) where i^sMS = ^'O- 

Since S, Sa and E are homogeneous functions of order 1 in iV, {dS/dE)]\f Y 
and its analogue {dSa/dE)^^ are homogeneous functions of order in iV so 
that 

l/{dS/dE)Ny = T + Oil/N), l/{dSaldE)N.y = T„ + 0(1/7V), (3) 

where T,Tq, represent the temperatures corresponding to S,Sa, respectively, as 
N ^ 00. The behavior in (3) is critical in understanding the entropy function 
slope at K, the Kauzmann point [S'sms = 0]. There arc two distinct possibilities. 
The slope at K is either finite, as shown explicitly in Fig. 1 and consistent with 
our rigorous analysis, or infinite. The former corresponds to a positive Tk, 
while the latter corresponds to an absence of a Kauzmann paradox. Almost all 
explicit calculations show the former as the usual behavior, two of which are 
presented below. 

We focus on the behavior of the growth function,^® which is adapted here 
with a slight modification. We divide the almost continuum energy spectrum 
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Figure 1: Schematic form of the generic entropy functions for various possible 
states. 
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E into bins, each of width Ae > 0, and index the bins by an integer j, with 
j = for the lowest energy bin, which contains all energies between Eq and 
So + Ae. We may take for Ae the average energy difference between any two 
consecutive energy levels of a single particle spectrum corresponding to its boimd 
state; its actual choice is not important for our discussion. Let W{B) denote 
the number of distinct configurations corresponding to energies lying in the bin 
with index B. The growth function-'^^ m(^) for fixed A^, V is given by JI^B) = 
W{B + 1)/W{B), and represents on an average the number of new configuration 
each older configuration generates as we add a quantum of energy to the system. 
It is easy to sec that /l(^) is a homogeneous function of order in A'^. Taking 
the logarithm on both sides and expanding S{B + 1) about B for finite A^, we 
immediately recognize that 

In 71(5) = (dS/dB) + 0{l/N) = Ae/T + 0{1/N) > (4) 

for non-negative T, where we have used (3). We immediately note that not only 
7Z(B) > 1, but also that it is monotone decreasing function of B. It achieves 
its maximum value JiiB) ^ oo at r=0, see point O in Fig. 1, and becomes 
finite as soon as T is raised above zero, decreasing monotonically to 1 at infinite 
temperatures. We similarly introduce the growth functions using W„, which 
satisfy Eq. (4) except that T is replaced by T^. To proceed further, we write 
W{B) = VFord(S)[l +^dis(S)/^ord(S)], where X« = WJW (0 < X„ < 1). It 
is easy to see that 

In /I = InTi^.d + + ^dis(/Idis/7*ord - !)]• (5) 

We emphasize that In /I and \njl^ take their limiting values Ae/T, and Ae/Ta, 
respectively, as N ^ oo; however, even for finite but large N, they are not 
far from their limits as they are homogeneous functions of order 0. For finite 
A'', Xdis is strictly non-zero but vanishes a,s N ^ oo. Let us apply (5) at 
B = Bk, which corresponds to the bin containing E = Ek > Eq. Since CR is 
the equilibrium state at E}^, we have, for A^ oo. In /J = In /lord = Ae/To^d, 
the CR growth function (at Ek), which must be finite since Ek corresponds to 
a non-zero CR temperature. The effect of /xjis noticed only for finite A'' as we 
discuss below. For finite A'^, In /I must be close to Ae/Tord- This immediately 
implies that /Idis cannot diverge at Ek ■ It must be finite, and must remain finite 
even as N ^ oo. For this to be so, it is clear that Tjis = 7k cannot be zero 
at Ek- Thus, as long as Ek > Eq, /idis must remain finite as A'^ ^ oo, and 
there must exist a positive Kauzmann temperature Tk- This constitutes the 
second rigorous proof of our central result. In CSM, 5 represents the classical 
configurational entropy.'' In QSM, S represents the total entropy." Thus, Tk 
corresponds to vanishing of different entropies in the two cases. Any attempt to 
estimate the classical configurational entropy in QSM, where it has no meaning, 
will require some sort of approximation, which we do not consider here as we 
are interested in rigorous results. Therefore, we will not pursue this point fur- 
ther. 
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Since the minimum energy state E = Ey^ in SMS corresponds to a positive 
temperature, the SMS state cannot exist below it. Hence, an ideal glass tran- 
sition must be invoked at this Kauzmann point (T = Tk). The discussion also 
establishes that the Kauzmann point in the disordered phase exists only when 
there exists another equilibrium state; otherwise, there will be no partitioning 
and, therefore, there will be no Kauzmann point; see also the explicit calculation 
below. 

We consider two CSM models in which we obtain positive Kauzmann tem- 
perature. The calculations are carried out exactly 

1. One-dimensional Polymer Model. The example also shows how p is 
used to distinguish different phases. Consider m-component spins Sj located 
at site i of a one-dimensional lattice with periodic boundary condition (Siv-i-i 
=Si ). Each spin can point along or against the axes (labeled 1 < fc < m) 
of an TO-dimensional spin space and is of length -^/m- The spins interact via a 
ferromagnetic nearest-neighbor interaction of strength (—J), with K = J/T > 
0. The high- temperature expansion of the PF is 

Zn {K, m) = ^K^m^, (6) 

where m > can be taken as a real number, and describes polymers with 

multiple bonds and loops, with K, and m denoting the bond and the loop 
activity, respectively.^'' The transfer matrix exp(iirS • S') has eigenvalues [x = 
exp{Km)] 

Xi=x + l/x + 2(m - 1), A2 = a; - 1/x, X3 = x + 1/x - 2, (7) 

that are 1-fold, m-fold, and (m — l)-fold, respectively.^* The eigenvalue Ai is 
dominant at high temperatures and describes the disordered phase. Its eigenvec- 
tor (^disl = (*l /\/2m, where (2fc| ( {2k + 1|) denotes the spin state pointing 
along the positive (negative) fc-th spin-axis, has the correct symmetry to give 
zero magnetization (p = 0). For m < 1, A2 becomes dominant at low tem- 
peratures [x > Xc = — m)] and describes the ordered phase, with a phase 
transition at Xc ■ The corresponding m eigenvectors are given by the combination 

(^ord = [(2^1 - {2k+l\]/V2, which are orthogonal to (^dis| ■ These eigenvec- 
tors have the symmetry to ensure p ^ 0. We note that both eigenvalues give 
the same free energy again at absolute zero. This means that if the eigenvalue 
Ai is taken to represent the metastable phase above a;c, its free energy must 
become equal to that of the equilibrium phase (given by A2) at absolute zero, in 
conformity with the proof given above. The bond and the loop densities 4>b = 
d\n\i/d\nK, and = d\nXi/d\nm are used to calculate the entropy per 
site s = InAi — <^b1ii^— (pLlnm. The proper convexity requirements {d(j)B/ 
d\nK > 0, dipL/ 9 In TO > 0) for (6) are always satisfied for Ai; see, for example, 
the behavior of (pB in Fig. 2, where we have taken m = 0.7. As the high- 
temperature disordered phase represents a physical system, it cannot give rise 
to a negative entropy; however, its metastable extension violates it as shown in 
Fig. 2, where its entropy becomes negative below Tk = 0.266, which is lower 
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than the transition temperature T^. As m 0, both Tk and (Tk < Tc) move 
down to zero simultaneously, and there is no Kauzmann point since there is no 
ordered state any more, as argued above. Thus, our exact calculation confirms 
our results. 

It should be noted that the eigenvalues Ai and A2 are independent of the size 
of the lattice. Therefore, they can be used to describe not only the disordered 
and ordered phases, but also the SMS, which is the continuation of the disor- 
dered phase, even for a finite N . Thermodynamic limit is not necessary. For 
finite N, Z^"' [K, m) < Z{K, m) and the inequality becomes an equality only 
as A'' — !■ 00 for the proper choice of a depending on the temperature. 

2. Binary Mixture Model. We now consider an anti- ferromagnetic Ising 
model (J > 0) in zero magnetic field on a Husimi cactus, which can be solved 
exactly. "'^^ The cactus can be thought as an approximation of a square lattice, 
so that the exact Husimi cactus solution can be thought of as an approximate 
solution of the square lattice model. The two spin states represent the particles 
of two species. There is a sublattice structure at low temperatures caused by 
the anti-ferromagnetic interaction; particles of one species are found on one of 
the two sublattices. We identify this ordered structure as a crystal. We allow 
two-, and three-spin (J' 7^ 0) interactions within the same square; the latter 
ensures that the melting transition is first order. The interaction energy is 

S=jY^SS' + J'^SS'S". (8) 

The first sum is over nearest-neighbor spin pairs and the second over neighboring 
spin triplets within each square. We set J=l to set the temperature scale. The 
model is solved recursively, as has been described elsewhere. We introduce 
partial PFs Zm{'\ or J,) depending on the states of the spin at the m-th cactus 
level. It represents the contribution of the part of the cactus above that level to 
the total PF. We then introduce the ratio Xm = Zm{'\)/[Zm{'\) + ■^m(i)])Which 
satisfies the recursion relation 

Xm = f{Xm,+l,Xm,+2,v)/[f{Xm+l,Xm+2,v) + f{ym+l,ym,+2,l/v)], where 

f{x,x',v) = x^x' /u'^v^ + 2xx'yv'^ + x^yv^ + u'^x'y'^ + 2xyy' + y^y' jv^ with 
u= e^^ ,v = e'^'^ , y = 1 — X and y' = 1 — x' . There are two kinds of fix-point 
solutions of the recursion relation, which describe the bulk behavior. 

In the 1-cycle solution, Xm = 1/2 is the same at all levels and exists at 
all temperatures. This solution corresponds to the disordered paramagnetic 
phase at high temperatures and the SMS below the melting transition. In the 
2-cycle solution, Xm alternates between two consecutive sites. This solution 
corresponds to the low temperature AF-ordered phase, which represents the 
CR at half occupation. For J' = 0.01, we find that Tm = 2.753, where there is 
a discontinuity in the entropy per site of 0.0168. The SMS below Tm represents 
SCL, whose entropy density, sec Fig. 2 vanishes at Tk = 1.132, and whose 
specific heat (not shown) remains positive with a maximum at T = 1.26. At 
absolute zero, the entropy per site Sscl — —0.3466. The CR and SCL free 
energies per site become identical (= 2J) at absolute zero in accordance with 
the claim above. Thus, the free energy diagram we obtain in this case is similar 
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Figure 2: The bond and the entropy densities. Both models show an entropy 
crisis at a positive temperature. 
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to that in Fig. 1. 

Note that the one-dimensional exact calculation is not a mean-field type 
calculation. Hence, our results are not of mean-field nature. We also note that 
there is no frustration in the Ising model, which shows the ideal glass transition. 
Both calculations confirm the general predictions. In particular, we conclude 
that if SMS exists, it gives rise to a positive Kauzmann temperature provided 
the ideal glass energy Ek > Eq, the ideal CR energy. 
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